The optimal cloning of quantum coherent states is non-Gaussian 
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We consider the optimal cloning of quantum coherent states with single-clone and joint fidelity as figures of 
merit. Both optimal fidelities are attained for phase space translation covariant doners. Remarkably, the joint 
fidelity is maximized by a Gaussian doner, whereas the single-clone fidelity can be enhanced by non-Gaussian 
operations: a symmetric non-Gaussian l-to-2 doner can achieve a single-clone fidelity of approximately 0.6826, 
perceivably higher than the optimal fidelity of 2/3 in a Gaussian setting. This optimal doner can be realized by 
means of an optical parametric amplifier supplemented with a particular source of non-Gaussian bimodal states. 
Finally, we show that the single-clone fidelity of the optimal 1-to-oo doner, corresponding to a measure-and- 
prepare scheme, cannot exceed 1/2. This value is achieved by a Gaussian scheme and cannot be surpassed even 
with supplemental bound entangled states. 

PACS numbers: 03.67.-a, 03.65.Ud, 42.50.Dv 



The no-cloning theorem states that there is no quantum 
apparatus capable of perfectly duplicating an arbitrary input 
state Ql. This is a direct consequence of the linearity of quan- 
tum mechanics and a fundamental difference between classi- 
cal and quantum information. This theorem enables one of 
the most promising applications of quantum information the- 
ory, namely secure quantum key distribution. Moreover, the 
impossibility of perfect cloning machines is intimately con- 
nected to other impossible tasks in quantum mechanics 01 . 

Soon after the observation of the no-cloning theorem as 
a fundamental feature of quantum mechanics the question 
arose how well an approximative cloning machine could 
work. For the case of universal cloning of finite dimen- 
sionalpure states this question was addressed and answered 
in|3H|l|i|l|i|l. There, the figure of merit was the fi- 
delity, i. e. the overlap between hypothetically perfect clones 
and the actual output of the imperfect doner. In particular, it 
was shown that judging single clones leads to the same opti- 
mal doner as when comparing the joint output with a tensor 
product of perfect clones ]6l|7j]. 

Recently, more and more attention has been devoted to con- 
tinuous variable systems, especially to states with Gaussian 
Wigner function - so called Gaussian states. Besides their out- 
standing importance in quantum optics and quantum commu- 
nication, in particular quantum cryptography fioll . they pro- 
vide a closed testbed within which many of the otherwise 
hardly tractable problems in quantum information become 
feasible. Restricting to the Gaussian world, i. e. to Gaussian 
operations on Gaussian states, led for instance to solutions to 
otherwise open problems in the theory of entanglement mea- 
sures Hill , quantum channels fl2ll and secret key distillation 
Jl3ll . Similarly, the problem of cloning in particular coherent 
states by Gaussian operations has been addressed in 
The obtained doner was shown to be optimal within the class 
of Gaussian operations by exploiting the connection with state 
estimation 1 16]. However, it remained unclear whether Gaus- 
sian operations really lead to the optimum, even under the as- 
sumptions typically made in the literature such as phase space 



translation covariance or output symmetry. 

The present Letter is concerned with the problem of opti- 
mally cloning coherent states without imposing any restric- 
tions on the cloning operation. After recalling some prelimi- 
naries, we will prove that w. 1. o. g. one can restrict to covari- 
ant doners, for which a powerful characterization will be pro- 
vided. Based on this, we will show that, in contrast to the finite 
dimensional case, the optimal doner depends on whether we 
judge single clones or test the clones jointly. Surprisingly, in 
the latter case the known Gaussian doners turn out to be op- 
timal, whereas with respect to the single-clone fidelity, non- 
Gaussian operations can perform better. 

The problem of finding the optimal doner reduces to find- 
ing the dominant eigenstate of an appropriate operator. For 
the optimal l-to-2 doner this eigenstate is directly linked to an 
optical implementation: it is the bimodal state of light that has 
to be injected on the idler mode of an optical parametric am- 
plifier and the input port of a beam splitter. We envision that 
a few-photon approximation of this doner, only sub-optimal 
but yet non-Gaussian, might be feasible, making it possible to 
experimentally demonstrate this fidelity enhancement. 

In addition, we will show that a l-to-00 doner based on a 
measure-and-prepare scheme cannot exceed a fidelity of 1/2, 
not even with supplemental bound entangled states. Extended 
discussions of the mathematical details 1 17] and the quantum 
optical aspects fiH will be reported elsewhere. 

Phase space and coherent states. Consider a system of 
n harmonic oscillators with respective canonical operators, 
or optical field quadratures, (Q\, Pi, . . . , Q n , P n ) =: R and 
the corresponding phase space 3 = M. 2n , which is equipped 
with an antilinear symplectic form cr(£, 77). Translations in 
this phase space are governed by the Weyl or displacement op- 
erators = e lcr ^' R \ £ <E S, which in turn obey the Weyl 
relations 
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implements the symplectic form via cr(£, 77) = ^ T -a-rj. Tensor 
products in Hilbert space correspond to direct sums in phase 
space, and in particular W%, = Wq j . where each of the 
£j G M 2 belongs to a single mode. 

The expectation values of the Weyl operators completely 
determine a state, and the resulting function, which is the 
Fourier transform of the Wigner function, is called the char- 
acteristic function. For a coherent state, it is a Gaussian of the 
form 

x(e)-tr[p^]- e -^ T -^- 8dr -«, (2) 

with covariance matrix 7=1 and displacement vector d. Co- 
herent states are translations of the harmonic oscillator ground 
state Wf|0) = |£) with d = a ■ £. In quantum optical set- 
tings, position and momentum coordinates correspond to the 
real and imaginary parts of the complex field amplitude. 

Figures of merit. The fidelity quantifies how close two 
states p\ and P2 are fl9ll . Here, we only consider the case of 
pure input states, so we can simply set f(pi, p-i) = tr[/?i pi\. 
A 1-to-n cloning transformation T (a "doner" for short) by 
definition takes systems in the pure input state p into n sys- 
tems whose state is close to n copies of p. We can express this 
by requiring the fidelity 

f joiIA (T,p)=tx[T(p)p® n ] (3) 

to be as large as possible. This is a very demanding criterion, 
as it also evaluates whether the clones are nearly independent. 
Instead, we might just evaluate the quality of an individual 
clone, say the i th , 

fi(T, p) = tr[T(p) (1 <g> ■ • • ® 1 ® p {i) <g> 1 ® ■ ■ ■ ® 1)] , (4) 

where the upper index denotes the position in the tensor prod- 
uct. Since a single such fidelity can be trivially put to one by 
copying the input onto the i' h clone, we have to maximize a 
weighted sum J^i Kfi(T, p) with positive weights A^. 

Further options arise from the choice of the set of states p 
that we want to clone optimally. Here, we consider the family 
of coherent states p = with |£) = W^|0). We define 

/joint(T') and fi(T) as the respective worst-case fidelities, i. e. 
the minima of (0 and @ over all coherent states p. Note 
that this is different from the usual case of universal doners 
in finite dimensional Hilbert spaces, where one considers the 
minimum with respect to all pure states. This is connected 
to the infinite number of dimensions of the continuous vari- 
able Hilbert space: Even minimizing (0 or (|4} over all pure 
squeezed Gaussian states (a larger though still very small sub- 
set of all states) would already yield a zero fidelity for all T. 

Our goal is thus to find the optimal worst-case joint fidelity 

/joint = SUp/joint(T) = SU P inf /joint (T, p) (5) 
T T p£ con 

as well as the convex set of achievable n-tuples of single-clone 
fidelities (/i(T), / 2 (T), . . . , f„(T)) as T varies overall don- 
ers. This is simplified by the fact that both fidelities are invari- 
ant under displacements in phase space, so we can choose the 



optimal doner to be covariant. Consequently, they are optimal 
with respect to both worst-case and average fidelities. 

Covariance. Let T be a 1-to-n cloning map. If displacing 
the input in phase space is equivalent to displacing the outputs 
by the same amount, then T is called (displacement) covari- 
ant: 

T(p) = Wf n] T {W i p W\) Wf n = T ( (p) (6) 

for all £ and p, where we have defined the shifted doner 
for later reference. The doners investigated in lfl4l FlUl were 
restricted to be covariant. However, this need not be assumed, 
but rather comes out as a property of the optimal doners. As 
in the case of cloning of finite-dimensional systems 0], the 
core of the argument is averaging over the symmetry group: 
we have, for / = / joint or / = £\ A*/i, respectively, 

f(T) = Mf(T, <M f /(T e ,|0)(0|) 

« (7) 
= /(M e T ?J |0)(0|) = /(M^). 

Here stands for "mean with respect to £" and is imple- 
mented by an invariant mean I20I1 . So, the averaged and thus 
covariant doner is at least as good as T for all T, and we can 
restrict the search to the covariant case. 

Note that the output of such doners could be singular for 
this phase space average. A detailed argumentation shows, 
however, that this is not optimal for the fidelities considered 
here 0. 

Optimizing covariant doners. In the Heisenberg picture, 
(the adjoint of) a covariant doner maps Weyl operators onto 
multiples of Weyl operators, 

T*(Wk,...t n )=t(S 1 ,...,{ n )W Et t i , (8) 

where £j is the pair of phase space variables of the i th clone. 
In terms of characteristic functions of input and output states, 
t acts as a characteristic function of the doner: 

Xoot(€l,-...fn)=*(6 ) ---,en)Xm(Ei^)- ( 9 ) 

The condition of complete positivity requires that t is the char- 
acteristic function of a state pr, plus a fixed linear transforma- 
tion liUl . We call a doner Gaussian if t has a Gaussian form 
and it thus maps Gaussian input states onto Gaussian output 
states. 

Since fidelities are linear in T, and hence linear in pr, they 
can be expressed as expectation values of linear operators: 

f{T,p)=tv\p T F). (10) 

The appropriate operators Fj i n t and Fi do not depend on T, 
which allows us to reduce the supremum of the left-hand side 
of J10I to finding the state px (hence the map T) correspond- 
ing to the largest eigenvalue of F. This is the core of our 
method. 

Physically, the state pt is directly related (up to a suitable 
symplectic transformation) to the bimodal state that needs to 
be injected on the idler mode of an optical parametric ampli- 
fier together with the input port of a beam splitter in order to 
realize the doner T (see below). 
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Optimal fidelities. Since by Eq.@ the maximum fidelities 
are reached by covariant doners, we can restrict the further 
discussion to a vacuum input state p = |0)(0|. For Gaus- 
sian input states, the operators F in Eq. dlOl are themselves 
Gaussian, so that the respective fidelities / are optimized by 
Gaussian pure states pr, hence by Gaussian doners T. Con- 
sequently, the joint fidelity / joint (T) = /joint (T, |0)(0|) = 
tr \pt i*]oint] is maximized by a Gaussian doner. The maxi- 
mum fidelity is given by the largest eigenvalue of an appropri- 
ately defined operator i*j int> that is 

sup /joint CO = max spec (Fj oint ) = -. (11) 
t n 

Thus, the unique optimal doner in this case is the know Gaus- 
sian doner of fllfllfl6ll. 

For the single-clone fidelity, we have to maximize the 

weighted sum J2i=i *i ft = tr iPT J27=i ^ F i\- Since a lin " 
ear combination of Gaussian operators does in general not 
have Gaussian eigenfunctions, it turns out that the optimal 
doners with respect to single-clone fidelities are not Gaussian. 
For simplicity, we restrict the following discussion to the 1-to- 
2 cloning problem. In this case the maximum of the weighted 
sum of single-copy fidelities Ai/i + A 2 / 2 = ti[pr F] is the 
largest eigenvalue of the operator 

F = Ai e -W?+i> 2 2 )/2 + A 2 e -(Q!+^ 2 )/2. (12) 

A simple numerical method to find this eigenvalue is to iter- 
ate 4> n +i = F(f> n /\\F(f) n \\. Varying the weights Ai yields the 
fidelity pairs (/1, / 2 ) along the solid curve in Fig.^ In com- 
parison, the best Gaussian doners are given by rotation in- 
variant Gaussian wave functions with appropriate squeezing, 
and the resulting fidelity pairs are plotted in Fig.^ as a dot- 
ted curve. At the intersection with the diagonal of symmetric 
fidelities lie the respective optimal doners. For the optimal 
non-Gaussian doner, we obtain /1 = /2 ~ 0.6826, which is 
strictly higher than the fidelity of the optimal Gaussian doner, 
namely h = h = 2/3 (cf. O). 

Studying doners which are described by highly squeezed 
non-Gaussian states pt reveals that on the curve of optimal 
fidelity pairs the points with fi = 1 and f 2 = 1 are ap- 
proached with infinite slope I17II . It is thus clear that the it- 
eration for the largest eigenvalue does not become singular. 
This regime is of potential interest in quantum key distribu- 
tion, since nearly perfect clones for the legitimate recipient 
combined with clones of non-trivial fidelity for the eavesdrop- 
per would be the hallmark of a successful cloning attack. On 
the other hand, the potential room for this regime is tiny as it 
is already proven that Gaussian attacks are optimal for a large 
class of quantum key distribution protocols where the channel 
is probed via second-order moments of the quadratures I22I1 . 

Optical implementation. The Gaussian symmetric doner 
can be realized by linear amplification of the input state, fol- 
lowed by distributing the output state into the two clones with 
a balanced beam splitter 12311 . This corresponds to the setup 
shown in Fig.|2]where the idler mode of the amplifier (bi) and 
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FIG. 1: Achievable pairs (/i,/a) of single-clone fidelities in l-to-2 
cloning of coherent states. The dots represent the optimal Gaussian 
doner, while the solid curve indicates optimal non-Gaussian opera- 
tions. Fidelities in the lower left quadrant are accessible to measure- 
and-prepare schemes. Classical mixtures of the two "trivial" doners 
fall onto the dashed line. The dash-dotted diagonal marks symmetric 
doners, with intersection points corresponding to the classical, best 
Gaussian, and optimal cloning, respectively. The inset shows the in- 
finite slope at /1 = 1 for non-Gaussian doners as opposed to the 
Gaussian case. 

the second input mode of the beam splitter (b 2 ) are both ini- 
tially in the vacuum state. Let us now analyze the cloning 
transformation that results from injecting an arbitrary two- 
mode state at modes b\ and b 2 . If the intensity gain of the 
optical parametric amplifier is 2, the modes where the two 
clones emerge are related to the input modes via the canonical 
transformation 

01 = am + Ol + b 2 )/V2 

I ' (13) 
a 2 = a in + {b\ - b 2 )/V2. 

From this expression, it is straightforward to check that the 
underlying doner is displacement covariant. Moreover, if the 
input is in the vacuum state p = |0) (0|, the single-clone fideli- 
ties amount to expectation values of the observables 

F = e -(Qi+C?2) 2 /4-(Pi-P2) 2 /4 

(14) 

F 2 = e -Wi-Q=) 2 /4-(Pi+P 2 ) 2 /4 

where (Qi,Pi) and (Q 2 ,P 2 ) are the canonically conjugate 
field quadratures of modes 6j and b 2 , respectively. This ex- 
actly coincides with expression H 2i up to a symplectic rota- 
tion, namely a beam splitter transformation. 

Consequently, the problem of finding the optimal doner 
reduces to finding the eigenstate with highest eigenvalue of 
X1F1 + X 2 F 2 , that is, to find the optimal bimodal state to 
be injected in modes b\ and b 2 . Note that if hM is an EPR 
state, i.e. a suitable infinitely squeezed state 12411 . then this 
corresponds to the two extreme points of the solid curve in 
Fig-E The symmetric case Ai = A 2 is obtained by choosing 

oc 

IV) = ^c„|2n)|2n) (15) 

n=0 
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FIG. 2: Optical scheme of a displacement-covariant doner. The in- 
put mode etj n is injected on the signal mode of an optical parametric 
amplifier (OPA) of gain 2, the idler mode being denoted as b\. After 
amplification, the signal mode is divided at a balanced beam splitter 
(BS), resulting in two clones in modes a\ and a^. The second input 
mode of the beam splitter is noted i>2- If both b\ and &2 are initially 
in the vacuum state, the corresponding doner is the Gaussian doner 

of QUI El- 

In contrast, if we inject a specific two-mode state 
|^) into b\ and 62, we can generate the whole set of displacement- 
covariant doners, in particular the non-Gaussian optimal one. 



where |n) are Fock states and the probability amplitudes c n 
correspond to the dominant eigenstate of F\ + F2. Trunca- 
tions of this state to finite photon numbers correspond to sub- 
optimal doners: keeping only the vacuum term ra = 0we get 
the optimal Gaussian doner with fidelities 2/3, while allow- 
ing for n < 2 yields the higher fidelities fi = /a ~ 0.6801 > 
2/3. The experimental realization of this doner does not seem 
unrealistic, given the recently proposed schemes for condi- 
tionally preparing arbitrary bimodal states of light based on 
linear optics l25Tl . In the limit n — ► 00, we arrive at the opti- 
mal doner with f\ = fa « 0.6826. 

Independent studies of the cloning fidelities of coherent 
states in finite-dimensional Hilbert spaces and their numeri- 
cal extrapolation has indicated that the optimal fidelity ranges 
between 2/3 and 0.699, which is in accordance with our result 

m. 

Optimal classical cloning. Let us finally consider a clas- 
sical l-to-00 cloning map T which is realized by measur- 
ing and repreparing the system. From the line of arguments 
above, T can be assumed to be covariant. Since composing 
this doner with time reversal r leads to a completely positive 
map, t o T*(W P:q ) = Xt(V2p, V%q) W- P , q with xt(p, q) 
the characteristic function of a state. Computing the fidelity 
for coherent input states immediately yields: 
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/clas S ical(r,|0}(0|) = -tr[p r |0}(0|] 
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,, < -■ (16) 
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The bound is reached by a heterodyne measurement and 
repreparation of coherent states, i. e. by a Gaussian scheme. 
This limit can not even be surpassed with the assistance of 
PPT bound entanglement I27I1 . since the respective maps are 
included in the above argumentation. In the case of an unas- 
sisted measure-and-prepare scheme an independent proof was 
recently given in 12811 . 
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